Example
> A = [0.96, 1.72; 2.28, 0.96];
> Asvd = svd(A)
sigma u vt
> Asvd.vt
matrix columns 1 thru 2
-0.8 -0.6
0.6 -0.8
> Asvd.u
matrix columns 1 thru 2
-0.6 -0.8
-0.8 0.6
> Asvd.sigma
vector elements 1 thru 2
3 1
> check = Asvd.u * diag(Asvd.sigma) * Asvd.vt
check =
matrix columns 1 thru 2
0.96 1.72
2.28 0.96
> a = rand(100,20);
> b = svd(a)
b =
sigma u vt
> whos(b)
Name Class Type Size NBytes
sigma num real 1 20 160
u num real 100 20 16000
vt num real 20 20 3200
Total MBytes = 0.019360
0
> c = svd(a,"A")
c =
sigma u vt
> whos(c)
Name Class Type Size NBytes
sigma num real 1 20 160
u num real 100 100 80000
vt num real 20 20 3200
Total MBytes = 0.083360
0
> d = svd(a,"N")
d =
sigma u vt
> whos(d)
Name Class Type Size NBytes
sigma num real 1 20 160
u num real 0 0 0
vt num real 0 0 0
Total MBytes = 0.000160
0